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PREFACE 


This  report  summarizes  the  numerical  work  performed  by  California 
Research  and  Technology,  Inc.,  under  subcontract  to  The  University  of  Dayton 
Research  Institute  on  Phase  II  of  USAF  Contract  No.  F33615-75-C-5052.  The 
.,  authors  are  Mr.  M.  Rosenblatt,  Mr.  G.  E.  Eggum,  Mr.  L.  A.  DeAngelo,  and 

Mr.  K.  N.  Kreyenhagen.  The  experimental  work  performed  by  The  University  of 
| Dayton  Research  Institute  on  Phase  IIs  entitled  "Impact  of  Soft  Bodies  on 

Jr  . / 

Jet  Engine  Fan  Blades",  is  presented  in  UDRI-TR-76-33  . and  will  be  published 
, as  an  AFML  Technical  Report.  Work  on  Phase  I of  this  contract  was  nresented 

in  AFML-TR-76-34 , entitled  "Ordnance  Impacts  on  Jet  Engine  Fan  Blades". 

The  effort  described  was  conducted  in  support  of  Project  No.  7351  and 
Task  735106  for  the  Metals  Behavior  Branch  of  the  Metals  and  Ceramics 
> Division  during  the  period  July  1975  to  December  1975.  The  contract  monitor 

was  Capt.  James  S.  Wilbeck. 


This  report  was  submitted  by  the  authors  in  August  1976  for  publication 
as  a Materials  Laboratory  Technical  Report. 
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1.1  PROBLEM  AREA  AND  BACKGROUND 


Foreign  objects  (for  example  rocks,  birds,  tools,  etc.) 
ingested  into  turbine  engines  can  cause  substantial  or  catas- 

% 

trophic  damage  to  the  blades.  This  problem  is  of  particular 
importance  with  respect  to  the  development  and  utilization 
of  new,  light-weight  blade  materials  which  may  be  quite 
vulnerable  to  damage  by  such  impacts. 

The  response  of  a turbine  blade  to  a foreign  object 
impact  can  be  separated  into  local  (or  short  time)  and  global 
(or  long  time)  response.  Local  response  involves  the  force 
and  impulse  loading  and  the  damage  processes  in  the  immediate 
vicinity  of  the  impact.  Global  response  involves  the  longer- 
term  response  of  the  entire  blade,  including  the  attachment 
to  the  hub. 

The  objective  of  the  exploratory  program  reported 
herein  has  been  to  investigate  the  local  impact  response  of 
blades  to  soft  body  impacts  (e.g.  birds)  using  numerical  and 
analytical  techniques.  Once  the  local  impact  force  and  impulse 
are  characterized,  these  can  be  used  as  loading  boundary  con- 
ditions for  separate  analyses  of  the  global  response  of  blades 
and  attachment  structures.  (In  another  part  of  this  program, 
bullet  impacts  on  titanium,  boron  aluminum,  and  graphite 
epoxy  plates  are  being  analyzed  in  both  the  normal  and  oblique 
impact  configurations.  That  work  will  be  described  in  a 
separate  report.) 


1.2  APPROACH 


Computational  codes  were  used  as  the  primary  tools 
to  investigate  local  soft  body  impact  phenomenology.  The 
codes  predict  the  time-resolved  pressure  and  velocity  fields 
^ in  the  impacting  object.  They  also  predict  the  spatial  and 

. temporal  details  of  the  impact  forces  acting  on  the  target, 

as  well  as  the  resulting  target  response.  Figure  1 shows  the 

i 

cases  treated;  these  were  chosen  to  qualitatively  and  quanti- 
tatively investigate  the  phenomenology  of  soft  body  impacts 
on  the  face  as  well  as  the  edge  of  turbine  blades.  In  all 
cases  the  impacting  object  was  a 3-oz  sphere  of  water 
(r  = 1.07  in.  or  2.73  cm)  with  a velocity  of  675  ft/sec.  In 
this  study,  water  was  used  to  simulate  the  impact  properties 
of  birds. 

Cases  A and  B treat  a normal  and  a shallow  angle  impact 

i 

on  a rigid  surface,  while  Cases  C and  D treat  the  corres- 
ponding impacts  on  a 0.136-in.  thick  stainless  steel  plate. 
Comparisons  of  these  cases  show  the  effects  of  plate  defor- 
, mation  on  the  impact  loading.  Cases  E and  F involve  a shallow 

angle  impact  on  the  leading  edge  of  a representative  turbine 
blade . 

Cases  A,  B,  C,  and  D were  analyzed  using  the  two- 
dimensional  WAVE-L  finite  difference  code.  WAVE-L  is  a 
HEMP-type1  Lagrangian  code.  The  initial  conditions  for  the 
; 25°  oblique  impact  cases  (Cases  B and  D)  were  treated  in  a 

moving  frame  of  reference  using  a Galilean  transformation; 
the  3-D  oblique  impact  configuration  is  thereby  simplified 
to  an  axisymmetric  (2-D)  geometry.  This  transformation 
and  the  rigid  and  SS403  plate  impact  cases  (A,  B,  C,  and  D) 
are  described  in  Sections  2.1  and  2.2.  Cases  E and  F were 
analyzed  using  the  three-dimensional  NONSAP  finite  element 
code2;  the  results  are  described  in  Section  2.3. 
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CASE  c CASE  D 


For  all  cases: 

Impact  velocity,  Vq  = 675  fps  (205.7  m/sec) 
Impacting  body  is  3 oz  water  sphere 

iJ  7 m = 85  gms , dQ  = 5.46  cm 

n 

CASE  E (Stainless  steel  SS403  blade) 

CASE  F (Boron-Aluminum  blade) 


Figure  1.  Conditions  for  Soft  Body  Impact  Analyses 


In  parallel  to  this  computational/analytical  effort, 

UDRI  has  performed  an  experimental  investigation  of  normal 
and  oblique  soft  body  impacts  on  rigid  and  thin  plates3.  The 
computational  and  experimental  results  are  compared  in 
Section  3. 

1.3  SUMMARY  OF  MAJOR  RESULTS 

Figure  2 illustrates  the  characteristic  behavior  of 
a water  sphere  impact  at  675  fps  onto  a rigid  or  semi-rigid 
surface.  Radial  flow  of  material  quickly  develops,  with 
radial  velocities  which  are  over  twice  the  impact  velocity 
occurring  near  the  plate  surface.  Figure  3 shows  pressure- 
radius  profiles  on  the  plate  surface  for  impact  Case  A at 
a sequence  of  times  after  impact.  Note  that  the  highest 
pressures  at  7 and  14  ysec  occur  near  the  outer  contact 
radius  between  the  flattening  water  sphere  and  the  plate 
surface . 

The  force  on  the  plate  is  the  integrated  area  under 
the  pressure- radius  profile.  Figure  4 shows  the  force  history 
on  the  plate  for  Case  A.  The  sharp  initial  rise  corresponds 
to  the  period  when  the  impact  shock  is  interacting  with  the 
droplet  flow.  It  is  followed  by  a relatively  long  period 
during  which  the  force  is  nearly  constant.  This  is  the 
period  when  most  of  the  axial  flow  diversion  to  radial  flow 
takes  place.  Finally,  the  force  begins  to  decay  as  the  rear 
of  the  drop  approaches  the  target.  (When  the  axial-radial 
turning  is  complete,  the  force  will  drop  to  zero.) 

Figure  4 uses  both  dimensional  and  non-dimensional 
(normalized)  coordinates  for  force  and  time.  The  time  normal- 
ization factor,  t,  is  the  time  required  for  the  sphere  to 
move  its  own  diameter  in  the  direction  perpendicular  to  the 
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Figure  2.  Velocity  Field  for  a Water  Sphere  Impact  at  675  fps  onto 
a Rigid  Surface  (Case  A)  at  41  and  157  ysec. 
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Figure  3.  Radial  Pressure  Profiles  on  Rigid  Surface  for 
Water  Sphere  Impact  at  675  fps  (Case  A). 
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Figure  4.  Force  History  Applied  to  Rigid  Surface  by  85  gm  Water  Sphere 
impact  at  675  fps  (Case  A) 


surface;  the  force  normalization  factor,  f,  is  the  average 
force  required  to  stop  the  momentum  perpendicular  to  the 
surface  in  time  t,  i.e. 


r 


M(VQsin0) ^ 
x D 

o 

where  1)Q  = initial  sphere  diameter 

VQ  = initial  velocity 
0 = obliquity  angle  (from  surface) 

M = sphere  mass 

Table  1 lists  these  factors  for  the  normal  and  oblique  impacts 
analyzed  in  this  study. 


V'osin0 

MV  sin0 
o 


TABLE  1.  NORMALIZING  FACTORS 


VQ  = impact  velocity  = 675  fps  = 20.6  x 10^  cm/sec 
Dc  = sphere  diameter  = 2.15  in.=  5.46  cm 


M = sphere 

mass  = 3 

o z = 8 5 gm 

Normal  Impact 
0 = 90° 

Oblique  Impact 
0 = 25° 

Velocity  normal 
to  surface 

VQsin0 

675  fps 

(20,600  cm/sec) 

285  fps 
(8,700  cm/sec) 

Momentum  normal 
to  surface 

MVqs in0 

1.74  gm-cm/psec 

0.735  gm-cm/psec 

Time  normalizing 

Do 

265  psec 

628  psec 

factor,  x 

Vs  in9 
o ~ 

Force  normalizing 

M(VQsin0) L 

14750  lbf 

2630  1 b f 
(1  . 17xl09  dynes) 

factor,  f 

Do 

(6 . 56x10^  dynes) 

Use  of  non-dimensional  (normalized)  time  and  force 
parameters  allows  the  results  of  the  numerical  solutions  to 
be  scaled  to  any  size  impacting  sphere,  and  perhaps  to  other 
impact  velocities  and  obliquities.  Figure  5 compares  the 
non-dimensional  force  histories  in  the  normal  and  25°  oblique 
impacts  onto  a rigid  surface  (Cases  A and  C) . The  relatively 
small  difference  between  the  two  cases  confirms  that  the 
normalizing  factors  are  reasonably  valid,  at  least  for  those 
impacts  where  the  velocity  component  perpendicular  to  the 
surface  is  between  285  and  675  fps. 

The  force  histories  applied  to  a rigid  surface  and 
to  a deformable  surface  (i.e.  a 0.136-in.  stainless  steel 
plate)  are  compared  in  Figure  6 (normal  impacts)  and  Figure  7 
(25°  oblique  impacts).  The  respective  rigid  and  deformable 
curves  have  the  same  general  shape  in  both  figures;  however, 
the  force  levels  are  lower  on  the  deformable  plates.  This 
is  because  the  plate  displaces  in  the  direction  of  the  sphere's 
velocity,  thereby  reducing  the  effective  impact  velocity  and 
the  subsequent  pressure  and  force. 

The  momentum,  or  impulse,  imparted  to  the  surface  is 
the  integrated  area  under  the  normal  force-time  profiles. 
(Friction  between  the  water  sphere  and  the  impacted  surface 
is  assumed  to  be  negligible,  so  impulse  is  imparted  only  perpen- 
dicular to  the  surface.)  Figure  8 shows  the  time-histories 
of  the  impulses  delivered  to  the  targets  for  Cases  A,  B,  C, 
and  D.  (Impulse  is  normalized  by  dividing  by  the  initial 
perpendicular  component  of  momentum  of  the  sphere,  MV'osin0). 

All  the  time  histories  are  fairly  similar  on  this  normalized 
plot,  with  the  deformable  plate  impacts  (Cases  C and  D) 
falling  somewhat  below  the  rigid  target  impacts  (A  and  B) . 
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Figure  5.  Comparison  of  Force  Histories  in  Normal  and  25°  Oblique  Impacts  of 
85  gm  Water  Spheres  Onto  Rigid  Surface  at  675  fps  (Cases  A and  B) 
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Steel  Plate  Due  to  675  fps  Impacts  of  85  gm  Water  Spheres  at  25°  Obliquity 

(Cases  B and  D) 


Perpendicular  Component  of  Normalized  Impulse  Delivered  to  Target 
versus  Time  for  675  fps  Impacts  of  85  gm  Water  Sphere 


For  the  cases  where  impacts  on  the  leading  edge  of 
blades  were  considered  (E  and  F) , a blade  segment  corresponding 
to  the  profile  of  the  J79  Stage  1 blade  at  70%  span  was  modeled 
using  the  3-d  NONSAP  finite  element  code.  Loading  on  the 
blade  was  derived  from  results  of  the  Case  B oblique  impact 
solution  (25°  oblique  impact  on  a stainless  steel  plate)  . 
Deflection  histories  were  monitored  at  several  stations  in 
the  blade  segment.  Figures  9 and  10  show  this  information 
for  selected  stations  in  stainless  steel  and  boron  aluminum 
blade  segments,  respectively. 
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Vertical  Displacement,  Az.  cm 


Vertical  Displacement,  Az,  cm 


SECTION  II 

THE  NUMERICAL  SOLUTIONS 


In  this  section,  the  numerical  solutions  defined  by 
Figure  1 in  Section  1.2  are  described  and  illustrated  by 
velocity  fields,  computational  grids,  pressure  contours,  and 
curves  showing  the  spatial  and  temporal  distribution  of 
forces  on  the  target.  The  properties  used  to  model  the  various 
materials  in  these  solutions  are  described  in  Appendix  A. 

2.1  IMPACTS  ON  A RIGID  SURFACE 

The  physical  phenomena  taking  place  during  the  impact 
of  a spherical  water  drop  on  a rigid  target  have  been  examined 
both  theoretically  and  experimentally  (for  example,  References 
4 and  5) . Although  there  is  not  a sharp  delineation  in  the 
physical  processes,  the  impact  phenomena  have  generally  been 
considered  to  occur  in  two  phases.  During  Phase  I of  the 
impact  event,  the  compressibility  of  the  liquid  drop  is  the 
dominant  characteristic.  High  pressures  occur  near  the  impact 
point  and  there  is  little  radial  flow.  Phase  II  is  dominated 
by  the  radial  flow  diversion  of  the  liquid  along  the  rigid 
boundary  at  velocities  which  can  be  over  twice  the  initial 
impact  velocity. 

For  normal  incidence  impacts,  the  nominal  peak  pressure 
which  occurs  during  Phase  I is  roughly  the  "water  hammer 


pressure : , l . e . , 


where 


Pw  = P cV 

M o n . 

pQ  = density  of  the  liquid  drop 

c = sound  speed  in  the  drop 

V = normal  component  of  the  impact  velocity 


For  the  normal  incidence,  675  fps  impacts  in  the  current 
study,  P^  = 2.8  kb. 


17 


Rarefaction  waves  from  the  droplet  free  surface  provide  relief 
from  the  water  hammer  pressure.  Because  the  impact  velocity 
in  the  current  study  is  low  compared  to  the  speed  of  sound 
in  water  (4500  fps) , the  pressures  during  Phase  I drop  very 
rapidly  (i.e.  within  a few  usee)  from  the  water  hammer  level. 

2.1.1  Case  A - Normal  Impact  on  Rigid  Plane 

Case  A treats  the  normal  impact  of  a water  sphere  on 
a rigid  surface  at  675  fps.  The  initial  grid  configuration 
for  this  numerical  solution  is  shown  in  Figure  11.  There  are 
20  computational  cells  across  the  sphere  radius. 

Figures  12  and  13  are  velocity  and  pressure  field  plots 
at  7 and  14  usee.  Figure  13  is  just  beginning  to  show  sig- 
nificant radial  deformation.  While  still  relatively  high 
(1-1.5  kb),  pressures  at  7 and  14  usee  have  already  dropped 
well  be'lew  the  water  hammer  pressure.  It  is  interesting  to 
note  that  the  peak  pressures  (Pmax>  1 kb)  within  the  drop  at 
this  point  do  not  occur  at  the  rigid  contact  surface,  but 
rather  in  a region  slightly  above  that  surface. 

Figures  14  through  18  show  the  development  of  high 
radial  velocity  flow  characteristics  (i.e.  Phase  II)  along 
the  rigid  surface.  The  peak  radial  velocity  is  about  twice 
the  impact  velocity.  The  location  of  the  maximum  pressure 
in  the  drop  moves  down  to  the  rigid  plane.  During  this 
period,  the  dynamics  are  essentially  governed  by  incompres- 
sible flow. 

Figure  3 (in  Section  1)  shows  pressure  profiles  along 
the  rigid  surface  at  several  times.  The  water  hammer 
pressure  (about  2.8  kb)  lasts  for  only  about  1 usee  for  the 
sphere  size  and  impact  velocity  involved  in  Case  A,  and  is 
thus  not  seen  in  Figure  3.  By  14  usee,  the  highest  pressure 
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Initial  Computational  Grid  for  85  gm  Water  Sphere  Impact  at  675  fpa 
on  Rigid  Surface  (Case  A) 
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Velocity  Field  and  Pressure  Field  at  7.14  ysec  After  85  gra  Water  Sphere  Impact 
at  675  fps  on  Rigid  Surface  (Case  A) 
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Velocity  Field  and  Pressure  Field  at  14.2  usee  After  85  gm  Water  Sphere  Impact 
at  675  fps  on  a Rigid  Surface  (Case  A) 


Velocity  Field  and  Pressure  Field  at  71.8  psec  After  85  gm  Water  Sphere  Impact 
at  675  fps  on  Rigid  Surface  (Case  A) 
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Figure  16.  Velocity  Field  and  Pressure  Field  at  102  psec  After  85  gm  Water  Sphere  Impact 
at  675  fps  on  Rigid  Surface  (Case  A) 


eld  at  129.5  usee  After  85  gm  Water 
on  Rigid  Surface  (Case  A) 
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is  about  1 kb.  At  both  7 and  14  psec,  the  peak  pressure 
occurs  near  the  periphery  of  the  contact  surface  between  the 
droplet  and  the  rigid  surface.  By  41  usee,  the  highest 
pressures  are  on  the  axis,  but  a local  maximum  occurs 
closer  to  the  periphery  of  the  contact  surface.  By  72  usee 
(which  corresponds  to  the  flow  field  in  Figure  15),  the 
pressure  distribution  on  the  rigid  plane  has  reached  a nearly 
uniform  distribution  out  to  about  2.5-cm  radius.  Beyond 
2.5-cm,  pressure  decays  rapidly  out  the  the  edge  of  the 
contact  area. 

Figure  19  shows  the  force  history  applied  to  tne  rigid 
plate  and  the  impulse  history  delivered  to  the  surface  for 
the  normal  incidence  impact.  (The  non-dimensional  time  and 
force  parameters  defined  in  Section  1.3  are  used  in  Figure  19 
and  in  subsequent  force  and  impulse  histories  in  this  report.) 

2.1.2  Case  B - 25°  Oblique  Impact  on  Rigid  Plane 

Case  B treats  a water  sphere  impacting  a rigid  surface 
at  25°  obliquity,  as  shown  in  Figure  20.  This  is  a 3-D  problem. 
However,  the  oblique  configuration  in  Figure  20  is  dynamically 
equivalent  to  the  configuration  in  Figure  21.  If  there  is 
no  friction  between  the  water  and  the  surface  (a  reasonable 
assumption) , the  problem  in  Figure  21  can  be  considered  as 
axisymmetric  insofar  as  the  impacting  sphere  is  concerned. 

The  sphere  will  apply  only  normal  forces  to  the  rigid  surface. 
Thus  a solution  to  the  impact  problem  of  Figure  21  can  be 
easily  manupulated  to  obtain  a solution  for  the  configuration 
in  Figure  20.  A Galilean  transformation  to  a moving  coordinate 
system  connects  the  two  configurations.  This  is  discussed  in 
Appendix  B. 
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Figure  19.  Normalized  Force  and  Impulse  Histories  on  Rigid  Surface  from  85  gm  Water 
Sphere  Impact  at  675  fps  (Case  A) 


V = 675  fps 
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Figure  20.  Geometry  of  675  fps  Impact  of  Water  Sphere  at  25°  Obliquity 
on  Stationary  Rigid  Surface  (Case  B) 


V = 675  sin  25 
= 285  fps 


Rigid  Surface 


V = 675  cos  25 
= 612  fps 


Figure  21.  Dynamically  Equivalent  Configuration  of  Normal- Incidence- 
Impact  of  Water  Sphere  on  Moving  Rigid  Surface  (Case  B) 


29 


The  phenomena  which  occur  during  Case  B are  fundamenta i ly 
the  same  as  those  which  occur  during  Case  A.  The  initial 
computational  grid  used  for  Case  B was  identical  to  that 
shown  in  Figure  11  for  Case  A.  Figures  22  to  25  show  the 
location  and  distortion  of  the  impacting  water  sphere  as 
well  as  the  pressure  fields  in  the  water  at  selected  times. 

Figure  26  shows  the  pressure  distribution  on  the  rigid  surface 
at  several  times,  and  also  the  contact  area  on  the  surface 
at  these  times.  Note  that  while  the  contact  area  is  still 
circular  in  this  oblique  impact  case,  the  center  of  the  contact 
area  moves  along  the  plate  with  a velocity  of  675  cos  25°=612  fps. 

Figure  27  shows  the  force  and  impulse  histories  for 

Case  B. 

2.2  IMPACTS  ON  A DEFORMABLE  PLATE 

The  conditions  used  for  the  impacts  on  a rigid  surface 
in  Cases  A and  B were  repeated  using  a 0.136- in  thick  stain- 
less steel  target  plate.  The  influence  of  the  deformation 
of  the  target  on  the  pressure  and  force  loading  functions  was 
of  primary  interest  in  these  calculations. 

2.2.1  Case  C - Normal  impact  on  0.136- in.  Stainless  Steel 
Plate 

Figure  28  shows  the  initial  conditions  for  this  solution. 
Figures  29,  30,  and  31  show  velocity  and  grid  deformation 
fields  in  the  water  drop  and  target  at  various  times.  The 
largest  plastic  deformation  (generalized  plastic  strain) 
experienced  by  the  end  of  the  solution  was  4.3  percent,  which 
occurred  on  the  rear  surface  of  the  plate  under  the  impact 
point  (i.c.  at  the  crown  of  tiie  bulge).  The  plate  will 
remain  deformed  after  the  droplet  interaction 
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Figure  22.  Pressure  Contours  at  49  psec  After  85  gm  Water  Sphere  Impact  at  675  fps  and 
25°  Obliquity  on  Rigid  Surface  (Case  B) 


Pressure  Contours  at  100  psec  After  85  gm  Water  Sphere  Impact  at  675  fps 
and  25°  Obliquity  on  Rigid  Surface  (Case  B) 
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Pressure  Contours  at  240  usee  After  85  gm  Water  Sphere  Impact  at  675  fps 
and  25°  Obliquity  of  Rigid  Surface  (Case  B) 


Figure  26.  Contact  Area  and  Pressure  Distribution  on  a Rigid  Plane  for  a 85  gm 
Water  Sphere  Impact  at  675  fps  and  25°  Obliquity  (Case  B) 
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Figure  27.  Normalized  Force  and  Impulse  Histories  on  Rigid  Surface  from 

85  gm  Water  Sphere  Impact  at  675  fps  and  25°  Obliquity  (Case  B) 
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Figure  23.  Initial  Computational  Grid  for  85  gm  Water  Sphere  Impact  at  675  fps 
on  0.136-in.  Stainless  Steel  Plate  (Cases  C and  D) 


i^ure  29.  Velocity  Field  at  65  psec  After  85  gm  Water  Sphere  Impact 
on  0.136-in.  Stainless  Steel  Plate  (Cases  C and  D) 
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Figure  30.  Computational  Grid  at  152  psec  After  85  gm  Water  Sphere  Impact  at  675  fps 
on  0.136-in.  Stainless  Steel  Plate  (Cases  C and  D) 
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Figure  31.  Velocity  Field  at  152  usee  After  85  gm  Water  Sphere  Impact  at  675  fps 
on  0.136-in.  Stainless  Steel  Plate  (Cases  C and  D) 


Figure  32  shows  stress  contours  in  the  drop  and  in  the  target 
at  64  psec.  Tensile  stresses  in  excess  of  2 kbar  (-29,000  psi) 
have  developed  on  the  rear  surface  of  the  target. 

Figure  33  compares  the  impulse  delivered  and  forces 
applied  to  the  target  versus  time  in  Cases  A (rigid  surface) 
and  C (0.136-in.  stainless  steel  plate).  It  is  seen  that  the 
deflection  of  the  plate  affects  the  magnitude  and  timing  of  the 
applied  forces,  with  the  forces  being  lower  in  the  deformable 
case . 

2.2.2  Case  D - 25°  Oblique  Impact  on  0.136-in.  Stainless 
Steel  Plate 

The  Case  B calculation  described  in  Section  2.1.2 
treated  a 25°  oblique  impact  on  a rigid  plate.  By  comparing 
Figures  26  and  3,  it  is  seen  that  pressure  profiles  in  the  25° 
oblique  impact  on  a rigid  surface  are  very  low  as  compared  to 
profiles  in  the  normal  impact  (i.e.  Case  B vs  Case  A).  This  is 
because  the  normal  component  of  the  velocity  for  the  25° 
obliquity  impact  is  only  285  fps.  Thus  it  was  expected  that 
the  plate  deflection  in  a 25°  impact  onto  a 0.136-in.  stain- 
less steel  ( deformable ) plate  would  be  much  smaller  than  seen 
in  the  previous  section  for  the  normal  impact  onto  the  same 
plate  (Case  C) . We  therefore  used  the  same  approach  in  Case  D 
to  analyze  the  25°  impact  on  a deformable  plate  as  had  been 
used  in  the  rigid  plate  oblique  impact,  (Case  B)  i.e.,  a 
Galilean  transformation  was  made  to  the  coordinate  system 
in  Figure  21,  in  which  the  sphere  strikes  the  plate  normally 
at  675  sin  25°  = 285  fps,  and  the  plate  moves  laterally  at 
675  cos  25°  = 612  fps.  The  motion  of  the  plate  is  ignored 
during  the  integration  because,  1)  shear  stresses  (friction) 
across  the  interface  are  assumed  to  be  negligible,  and  2)  the 
deflection  of  the  plate  is  assumed  to  be  small. 
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Using  these  assumptions,  the  velocity,  grid,  and 
pressure  fields  shown  in  Figures  34  to  37  for  various  times 
were  obtained.  The  deflection  of  the  plate  is  indeed  seen 
to  be  small,  i.e.,  the  maximum  deflection  is  0.162  cm  or  3%  of 
the  droplet  diameter . This  supports  the  simplifying  assump- 
tions which  were  used  during  the  numerical  integration. 

Figure  38  compares  the  impulse  force  histories  applied  by 
oblique  impacts  to  the  rigid  surfaces  (Case  B)  and  the 
deformable  surface  (Case  D) . Again,  the  forces  on  the  deform- 
able plate  are  somewhat  smaller  than  on  the  rigid  plate. 

2.3  IMPACTS  ON  THE  LEADING  EDGE  OF  TURBINE  BLADES 

For  these  analyses,  a 6- in.  segment  of  the  J79  Stage  1 
Turbine  blade  profile  at  70%  span  was  modeled  using  3-D 
NONSAP2.  The  blade  segment  was  rigidly  attached  at  both  ends, 
and  the  impact  occurred  on  the  leading  edge  at  the  center  of 
the  segment,  as  indicated  in  Figure  39.  Cases  E and  F are 
identical  except  for  the  turbine  blade  material  (Stainless 
Steel  in  Case  E and  Boron  Aluminum  in  Case  F) . 

In  the  NONSAP  model,  the  blade  was  uniformily  par- 
titioned in  the  x and  y-direction  into  6x6  elements,  with 
one  element  defined  through  the  thickness.  Each  element  is 
defined  by  12  nodes  as  illustrated  in  the  following  sketch: 
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CALIFORNIA  RESEARCH  AND  TECHNOLOGY  HAVE-1 
67S  FPS  IMPACT  ON  SSUC3  PLATE  AT  25  DEG  OBLIOUITT 


rurc  34.  Velocity  Field  at  68  |isec  Alter  85  gm  Water  Sphere  Impact  at  675  Cps  and 
25u  Obliquity  on  0.136-in.  Stainless  Steel  Plate  (Case  D) 
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Figure  36.  Velocity  Field  at  144  psec  After  85  gm  Water  Sphere  Impact  at  675  fps  and 
25°  Obliquity  on  0.136-in.  Stainless  Steel  Plate  (Case  D) 
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Figure  37.  Stress  Field  at  68  usee  After  85  gm  Water  Sphere  Impact  at  675  fps 
and  25°  Obliquity  on  0.136- in.  Stainless  Steel  Plate  (Case  D) 


For  these  solutions,  the  loading  from  the  bird  impact 
was  represented  by  a pressure  function  derived  from  the  2-D 
finite-difference  solution  described  in  Section  2.1.2  of  a 
water  sphere  impacting  at  67S  fps  on  a rigid  surface  at  25° 
obliquity  (Case  B) . Essentially,  it  was  assumed  that  the 
leading  edge  neatly  cuts  the  sphere  in  half  along  tne 
diametral  plane  which  is  shown  as  a dashed  line  in  Figure  39. 
The  blade  is  impacted  only  by  the  upper  hemisphere.  The  impact 
load  from  this  hemisphere  is  determined  from  the  pressure- 
radius-time  curves  from  Case  B (e.g.  in  Figure  26),  taken 
around  a 180°  arc.  The  contact  surface,  and  hence  the  impact 
load,  sweeps  across  the  top  surface  of  the  blade  at  Vocos0  = 
675  cos25°  = 612  fps 

The  pressure  loading  was  input  to  the  NONSAP  program 
as  forces  on  the  nodes.  The  forces  on  each  of  the  nodes 
were  calculated  at  a given  time  by  applying  the  appropriate 
loading  pressure  at  the  node  location  to  the  area  associated 
with  the  node.  The  resulting  force  was  applied  perpendicular 
to  the  blade  surface  at  the  location  of  the  node.  Thus,  the 
forces  from  the  Case  B solution  were  rotated  to  account,  at 
least  to  first  order,  for  the  fact  that  the  blade  is  not  flat. 

The  blade  materials  in  Cases  E and  F were  characterized 
as  linear  elastic  solids.  The  values  defining  the  material 
properties  are  listed  below: 

Young's  Poisson's 

Case  No . Material  Density  Modulus  Modulus 

E SS  403  7.86  gm/cm3  1.965  Mb  0.285 

F Boron-Alum.  2.40  gm/cm3  1.379  Mb  0.31 
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2.3.1  Stainless  Steel  (SS  403)  Blade 

Figures  40  to  43  show  vertical  (z -d i rect ion)  displace- 
ment as  a function  of  time  at  various  stations  on  the  stain- 
less steel  blade.  Negative  deflections  indicate  that  the 
station  is  depressed,  positive  deflections  mean  that  the 
* station  has  risen.  The  following  sketch  shows  the  locations 

of  output  stations. 

| 


The  maximum  deflection  for  the  SS403  blade  was  about 
0.2  cm,  which  occurs  at  the  impact  point  at  about  2S0  usee. 
Figure  40  shows  4 stations  along  the  center  of  the  blade. 
Stations  1,  2,  and  3,  which  are  near  the  impact  point,  have 
initial  downward  displacements.  Station  4,  on  the  trailing 
edge  of  the  blade,  first  displaces  upward,  reaching  an  up- 
ward maximum  displacement  at  roughly  200  usee. 

Figure  41  shows  displacement  histories  of  various 
stations  along  the  leading  edge.  Initial  displacement  is 
always  downward,  and  the  time  to  peak  deflection  is  consistently 
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Figure  40.  Vertical  Displacement  Histories  at  Selected  Points 
Downstream  of  Impact  of  85-gm  Water  Sphere  at  25° 
Obliquity  on  Leading  Edge  of  Stainless  Steel  (SS403) 
Blade  (Case  E) 
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Figure  43.  Vertical  Displacement  Histories  at  Selected  Points 
on  Surface  of  Stainless  Steel  (SS403)  Blade  After 
Impact  of  85-gm  Water  Sphere  at  25°  Obliquity  on  Leading 
Edge  (Case  E) 
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250  ysec  for  these  stations.  Figure  42  shows  stations  on  the 
rear  edge.  The  initial  displacements  are  upward,  with  a 
maximum  upward  deflection  of  about  0.04  cm,  occurring  in 
less  than  200  ysec.  This  upward  motion  along  the  trailing  edge 
is  smaller  than  the  subsequent  peak  downward  motions.  Figure 
43  shows  several  interior  stations.  The  stations  further 
from  the  impact  point  have  smaller  peak  deflections.  Also 
note  that  the  period  seems  to  increase  at  points  further  from 
the  impact. 

2.3.2  Boron-Aluminum  Blade 

Figures  44  to  47  show  displacement - 1 ime  histories  in 
the  Boron-Aluminum  blade  segment.  Quantitatively,  the 
results  are  similar  to  the  stainless  steel  blade;  however, 
the  peak  deflections  are  somewhat  larger.  At  the  impact 
point,  the  Bo-Al  peak  deflection  is  0.3  cm  as  compared  to 
0.2  cm  for  stainless  steel.  Also,  the  timing  is  somewhat 
different.  The  times  to  first  maximum  deflection  are  smaller 
by  roughly  50  usee  on  both  the  leading  and  trailing  edges 
of  the  Bo-Al  blade. 
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Figure  45.  Vertical  Displacement  Histories  at  Selected  Points 
Along  Leading  Edge  of  Bo-Al  Turbine  Blade  Impacted 
on  Leading  Edge  by  85-gm  Water  Sphere  at  25°  Obliquity 
(Case  F) 
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Figure  46.  Vertical  Displacement  Histories  at  Selected  Points 
Along  Trailing  Edge  of  Bo-Al  Turbine  Blade  Impacted 
on  Leading  Edge  by  85-gm  Water  Sphere  at  25°  Obliquity 
(Case.  F) 
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Figure  47.  Vertical  Displacement  Hisotires  at  Selected  Points 


on  Surface  of  Bo-Al  Turbine  Blade  Impacted  on  Leading 
Edge  by  85-gm  Water  Sphere  at  25°  Obliquity  (Case  F) 


SECTION  III 

COMPARISONS  WITH  EXPERIMENTS 

UDRI  has  conducted  an  experimental  investigation  of 
soft  body  impacts  on  nominally  rigid  and  deformable  targets 
at  normal  and  oblique  incidence3.  Recent  investigations6 
have  involved  transducer  measurements  of  pressure- t ime  his- 
tories at  stations  near  the  impact  point  in  the  target  plate. 
The  measured  histories  will  be  compared  with  corresponding 
calculations  in  this  section. 

Figure  48  summarizes  the  impact  conditions  for  the 
experiments  and  calculations.  Although  these  conditions  are 
nominally  similar,  significant  differences  in  impact  velocity, 
projectile  material  and  density,  and  projectile  shape  are 
apparent . 

3.1  NORMAL  IMPACTS 

Figure  49  compares  pressure  histories  at  the  impact 
point.  The  calculated  pulse  rises  very  steeply  to  about  the 
expected  water  hammer  pulse  (P^=  P0VQC  = 2.8  kb),  then  decays 
rapidly.  The  measured  pulse  rises  more  slowly,  attaining 
a maximum  after  40  Msec,  then  decaying  fairly  rapidly. 

(The  exact  maximum  was  not  determined;  it  was  stated  to  be 
in  excess  of  1.5  kb.)  After  about  50  usee,  both  the  calcu- 
lated and  measured  pressures  decay  more  gradually,  eventually 
approaching  the  dynamic  pressure  level,  P = 1/2  pVQ2~  0.2  kb. 

The  initial  differences  between  the  calculated  and 
measured  waveforms  may  be  due  to  the  limited  frequency 
response  of  the  transducers  (100  kHz)  , which  would  preclude 
resolution  of  the  initial  water  hammer  pulse  or  subsequent 
rapid  decay. 

The  oscillations  in  the  numerical  results  can  be 
smoothed  if  more  artificial  viscosity  is  used  in  these 
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Figure  48.  Summary  of  Impact  Conditions  for  Calculations  and  Experiments 


Figure  49.  Comparison  of  Calculated  and  Measured  Pressure  Histories  at 
Impact  Point:  Normal- Incidence  Impacts  on  Rigid  Surface 


type  analyses.  Only  a minimum  of  such  viscosity  was  used  in 
Case  A,  so  as  to  allow  the  early  pressure  spike  to  more 
accurately  develop. 

Figures  50  and  51  compare  measured  and  calculated 
pressure  histories  at  stations  1/2-in.  and  1-in.  from  the 
impact  point,  respectively.  The  experimental  data  lias  been 
shifted  in  time  to  agree  with  the  expected  time  delay  assoc- 
iated with  stations  at  1/2-in.  and  1-in.  from  the  impact. 

The  comparisons  for  these  stations  are  very  good. 

3.2  OBLIQUE  IMPACTS 

The  oblique  impact  comparisons  between  calculated 
and  measured  pressure  pulses  are  not  as  meaningful  because 
of  the  distinct  differences  in  projectile  shape  (sphere  vs 
right  cylinders)  and  the  velocity  variation  between  the  two 
tests.  Figure  52  shows  the  pressure  pulses  measured  at 
three  stations:  at  the  impact  center  (i.e.,  intersection  of 

the  cylinder  axis  with  the  target  surface) , and  at  points 
1-in.  "upstream"  and  "downstream"  from  that  point.  None  of 
the  pulses  show  a water  hammer  spike.  All  rise  relatively 
slowly,  and  remain  fairly  constant  thereafter.  However, 
the  pressure  levels  are  substantially  different.  The  up- 
stream station  attains  a pressure  level  of  about  110  bars. 


This  is  below  the  dynamic  pressure  corresponding  to  the  full 
velocity  for  the  gelatin  impact,  i.e.,  P = l/2pVQ2  = 180  bars), 
but  it  is  far  above  the  dynamic  pressure  corresponding  to 
the  normal  component  of  the  velocity,  P = l/2p  (VQsin0)2  = 32 
bars.  The  center-line  and  downstream  stations  attain 
pressures  which  correspond  closely  to  the  latter,  i.e., 

P = l/2p  (VQsin0) 2 . This  behavior  is  qualitatively  consistent 
with  the  tentative  observation  in  Reference  7 that  dynamic 


ure  transducer  data,  UDRI  Test  5380  (630  fps  gel 
CRT  calculation  of  Case  A (675  fps  water) 


lculated  and  Measured  Pressure  Histories  in  Target 
(1.5  in)  from  Impact  Point:  Normal-Incidence 


Figure  51.  Comparison  of  Calculated  and  Measured  Pressure  Histories  in  Target 
Surface  25.4  cm  (1  in)  from  Impact  Point:  Normal-Incidence 
Impacts  on  Rigid  Surface 


pressures  in  substantially  oblique  impacts  of  cylinders 
approach  l/2pVQ2  at  points  on  the  upstream  side  of  the  stag- 
nation point  (which  is  upstream  of  the  center-line  intersection), 
and  l/2(Vosin0)2  on  the  downstream  side. 

Figure  53  compares  the  calculated  pulse  (for  a spher- 
ical water  drop  impact-Case  B)  with  the  measured  pulse  (for 
the  gelatin  cylinder)  at  the  upstream  station.  The  waveforms 
are  qualitatively  very  different.  The  calculated  pulse  rises 
sharply  to  about  400  kbars,  then  decays  gradually,  reaching  the 
20-40  bar  level  by  15  ysec.  The  measured  pulse  climbs  grad- 
ually to  100-110  bars,  and  remains  nearly  constant  at  that 
level  for  300  ysec.  These  different  behaviors  indicate  that 
geometry  of  the  impacting  body  has  a potentially  important 
influence  on  the  loading  pulse. 
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Figure  53.  Comparison  of  Calculated  and  Measured  Pressure  Histories 
in  Target  Surface:  25°  Oblique  Impacts  on  Rigid  Surface 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  quantitative  and  qualitative  results  of  this  study 
demonstrate  the  detailed  nature  of  the  temporal  and  spatial 
descriptions  which  can  be  obtained  with  numerical  codes  for 
the  dynamic  impact  phenomena  involved  in  soft  body  impacts. 
Normal  impacts  on  rigid  or  deformable  plates  can  be  treated 
from  first  principals  using  the  2-D  finite  difference  WAVE-L 
code.  Oblique  impacts  are  more  difficult  to  analyze  because 
these  impacts  are  generally  3-D  in  nature.  For  the  case  of  a 
nominal ly- spherical  body  obliquely  impacting  an  essentially 
frictionless  surface  which  does  not  substantially  distort 
during  loading,  a Galilean  coordinate  transformation  (see 
Appendix  B and  Section  2.2)  can  be  used  to  transform  the  3-D 
oblique  impact  into  an  approximately  equivalent  2-D  axisym- 
metric  impact  problem.  This  transformation  technique  was 
applied  to  oblique  impacts  of  spherical  projectiles  into 
rigid  and  semi-rigid  targets.  For  oblique  impact  cases  in- 
volving large  target  deformations,  the  technique  is  not 
applicable . 

The  time-resolved  force  loading  functions  from  675  fps 
normal  and  oblique  soft  body  impacts  on  rigid  and  0.136-in. 
stainless  steel  plates  can  be  used  to  load  complete 
3-D  finite  element  models  of  turbine  blades,  including 
the  root  design.  Cases  E and  F in  this  study  examined  the 
local  behavior  in  3-D  of  a section  of  a turbine  blade  for  an 
edge-on  impact  at  25°.  A similar  approach  could  be  taken  using 
a more  complete  3-D  finite  element  blade  description. 

Potentially  important  differences  were  seen  between 
calculated  loading  pulses  for  obi iquely- impact ing  spheres, 
and  measured  pulses  for  obi iquely- impact ing  cylinders.  The 
differences  in  these  relatively  low-velocity  impacts  are  thought 
to  be  due  to  projectile  geometry  effects. 
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The  following  recommendations  are  made  for  additional 
efforts  directed  toward  understanding  the  soft  body  (bird) 
impact  problems: 

1 . Develop  realistic  material  model s for  birds  - 

The  numerical  codes  can  be  used  to  develop  and/or 
evaluate  material  and  loading  models  for  the  soft  body  impacts 
on  aircraft  components.  In  particular,  the  porosity  muscle/ 
skin  strength,  and  bone  structure  of  birds  may  influence  impact 
loading  and  local  damage  mechanisms.  Code  solutions  provide 
a precise  and  controllable  approach  for  understanding  the 
role  of  these  properties,  and  can  thereby  lead  to  development 
of  simpler  models. 

2 . Investigate  the  importance  of  projectile  shape  factors  - 

Codes  can  directly  evaluate  differences  in  the  force 
loading  on  a target  for  various  shapes  corresponding  to 
possible  bird  impacts.  This  information  is  needed  not  only 
to  evaluate  the  range  of  impact  loading  which  can  be  expected 
for  different  birds  and  impact  orientations,  but  also  to  help 
design  realistic  simulants. 

3 . Develop  a rapid  technique  for  using  the  early  impact 
phenomena  and  force  loading  descriptions  - 

This  is  needed  to  allow  for  efficient  3-D  finite 
element  code  prediction  of  survivability  of  turbine  blade 
designs . 
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APPENDIX  A 


i 


i 


MATERIAL  PROPERTIES 


1.  MATERIAL  MODEL  FOR  TYPE  403  STAINLESS  STEEL  (SS403) 

The  material  properties  for  SS403  described  below  have 
been  chosen  to  match  as  closely  as  possible  the  properties  of 
the  steel  being  used  in  the  experimental  work3,  which  was 
specified  by  the  supplier  to  have  an  ultimate  tensile  strength 
of  75520  psi,  and  a yield  stress  of  36540  spi. 

The  stress  levels  encountered  in  the  present  effort  are 
such  that  an  elastic-plastic  material  model  is  appropriate. 
Within  the  yield  surface,  deformations  are  considered  to  be 
elastic;  excursions  on  the  yield  surface  are  made  in  accor- 
dance with  the  Prandtl- Reuss  flow  rule.  The  location  of  the 
yield  surface  for  a material  element  is  a function  of  the 
local  generalized  plastic  strain  0,  e , where 


J\ J3  1J  1J 


and  e?j  is  the  plastic  strain  rate.  Work  hardening  is  thus 
included  in  the  model.  The  elastic  constants  used  in  the  model 
are  given  in  Table  A-l. 

TABLE  A-l.  ELASTIC  CONSTANTS  FOR  TYPE  403 
STAINLESS  STEEL 

Young's  Modulus 
Poisson's  Modulus 
Shear  Modulus 
Bulk  Modulus 

l 


28.5  x 

106 

psi 

0.285 

11.1  X 

106 

psi 

22.1  x 

106 

psi 

Figure  A- 1 shows  a simple  tension  stress -strain  curve 
for  stainless  steel  taken  from  Reference  8.  The  yield  stress 
(corresponding  to  0.2%  offset)  is  close  to  the  value  given 
for  the  steel  used  in  the  experiment.  This  curve  was  monoton- 
ically  extended  to  the  point  of  failure,  which  is  taken  to  be 
75520  psi  at  30%  strain.  For  efficient  use  in  the  solutions, 


the  curve 

has  been  fit  by 

a line 

segment  and  two  parabolic 

l 

segments , 

as  shown  in  Figure  A-2. 

The 

equations  for  the 

I 

analytical 

fit  are: 

a 

= aQe  for 

0<  e < 0.9875 

a 

* -V  bi2  - 

4aAci 

-e) 

/2a1for  0.9875<e<  .33 

0 

■ -b2*  b22  - 

4a2(c2 

-e) 

/ 2a?for  0 . 33  < e<  .30 

» 

where  aQ=2 

.857  x 1 03  psi 

al  = 1 

.472  x 10'9  psi" 

? 

, a2  = l 

.871  x 

-8  . -2 
10  psi 

bl  = 7 

.012  x 10  3 psi 

] , b2  = 

-1.335 

x 10  3 * * psi  ^ 

1 

cl  = 

. 9202  , 

C 2~ 

24.14 

2.  WATER  EQUATION  OF  STATE 

The  numerical  fit  of  Walker  and  Sternberg9  is  used  for 

the  equation  of  state  of  water.  This  equation  of  state  of 

water  is  valid  to  about  250  kbars.  Figure  A-3  shows  the 

Water  Hugoniot  up  to  50  kbars. 
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Figure  A-2.  SS403  simple  tension  stress-strain  curve  fit  used 
in  the  numerical  solutions  showing  the  comparison 
with  the  base  curve  from  Watter  and  Lincoln  (Ref.  8) 
and  the  extrapolation  to  large  strains 


76 


Pressure,  kb 


APPENDIX  B 


GALILEAN  COORDINATE  TRANSFORMATION 

A physical  event,  such  as  the  impact  of  a soft  body 
on  a plate,  cannot  be  influenced  by  the  particular  coordinate 
frame  used  by  an  observer  of  the  event.  The  usual  coordinate 
frame  of  reference  may  be  called  "the  laboratory  frame  of 
reference",  i.e.,  the  event  is  viewed  from  a stationary  frame 
of  reference  with  respect  to  the  earth.  A frame  of  reference 
(or  coordinate  system)  moving  with  a constant  velocity  with 
respect  to  the  laboratory  frame  of  reference  can  be  used  just 
as  well  to  describe  the  physical  event.  The  relationship 
between  a laboratory  and  moving  frame  of  reference  is  illus- 
trated in  the  sketch  below. 


For  simplicity,  and  for  the  application  in  Section  2.1,  we 
have  assumed  that  the  coordinate  axes  in  the  laboratory  and 
moving  frames  of  reference  are  parallel  with  the  relative 
motion  in  the  x direction;  that  is,  the  primed  system  is 
moving  in  the  x direction  with  speed  v relative  to  the 
unprimed  system. 


Now  if  the  two  origins  0 and  O'  are  coincident  at  time 
t = 0,  then  the  distance  00'  is  equal  to  vt . The  equations  of 
transformation,  according  to  classical  or  Newtonian  kinematics, 
are  then 

x = x ' + vt 

y = y' 

Z = Z ' 

t = t' 

where  x,  y,  z,  t,  are  in  the  laboratory  system,  and  x',  y', 
z',  and  t'  are  the  coordinates  of  the  moving  frame  of  reference. 

Velocities  in  the  two  frames  of  reference  are  then 
related  by 

x = x ' + v 

y = y' 

z = z ' 

Thus  the  dynamic  results  of  a physical  event  in  one 
coordinate  frame  is  sufficient  to  completely  determine  the 
dynamics  in  the  other  frame  of  reference. 
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